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ABSTRACT 



We use recently compiled position and velocity data for the globular cluster and planetary nebula 
subsystems in NGC 5128, the nearby giant elliptical, to search for evidence of past dwarf-satellite 
accretion events. Beyond a 10' (~ 11 kpc) radius in galactocentric distance, we find tentative evi- 
dence for 4 subgroups of globular clusters and 4 subgroups of planetary nebulae. These each have 
> 4 members within a search radius of 2' and internal velocity dispersion of < 40 km s~ x , typical 
parameters for a dwarf galaxy. In addition, 2 of the globular cluster groupings overlap with 2 of the 
planetary nebulae groupings, and 2 subgroupings also appear to overlap with previously known arc 
and shell features in the halo light. Simulation tests of our procedure indicate that the probability of 
finding false groups due to chance is < 1%. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: individual (NGC 5128) — galax- 
ies: interactions — galaxies: star clusters — globular clusters:general — planetary 
nebulae : general 



1. INTRODUCTION 

In the last decade or so, the presence of stellar 
streams within nearby galaxies has provided direct ev- 
idence for the continual build-up of massive galax- 
ies from the accretion of small neighbouring satellites. 
For example, in the Milky Way, the r e is the Sagit- 
tarius stream (llbata et al.1 119971 l2001ah iMaiewski etaTl 
120031: iMartmez-Delgado et all 120041 ). and in M31, a 
number of new streams have been recently discov- 
ered within the Pan- And romeda Archelogical Survey 
(|McConnachie et all I2009D . adding to the list of a l- 
ready known streams in that galaxy (jlbata et al. 2001b). 
As well, striking surface brightness features tracing a 
galaxy's merger hi story have b een found in NGC 590 7 
(IShang et al.1 11998ft. NGC 891 (IMouhcine et al.l [2010ai). 
NGC 4093 (IMartmez-Delgado et all 120081) . NGC 253 
and NGC 52 36 (IMalin k HadlevHl997l) . and UGC 10214 
(|Forbes et alll2003D . among others. Past work has shown 
that it is possible to use streams to trace the orbits of 
satellite galaxies that had merged with the Milky Way 
in the distant pa st (iLvnden-Bell k Lvnden-Belll Il995j ; 
lYoon fc Led [2001 . ILvnden-Bell k Lvnden-Bel (|l995T ) 
suggest that the tidal debris stripped out of the parent 
dwarf galaxy will not undergo much dynamical friction 
because its mass content will be low. 

In systems beyond our Local Group galaxies, finding 
these traces of old satellites becomes more difficult. But 
in addition to the field-star integrated light, it is pos- 
sible to use other old objects such as globular clusters 
(GCs) and planetary nebulae (PNe) to provide evidence 
for the presen ce of accretion remnant s dForte et alJll982t 
lMuzziol ll987l: iCote et alllliM 12001 IHilker et al.lll999f) . 
As well, iPipino et all (|2007ft have modelled the GC as- 
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sembly history in massive galaxies, and showed that in 
addition to the GCs that formed along with the galaxy 
halo, a significant fraction of GCs that make up the mas- 
sive galaxies were likely accreted from dwarf galaxies in 
later times. Indeed, subgroups of GCs have been found 
within the halos of massive galaxies or have been associ- 
ated with the stel l ar str e ams in the Sagittariu s stream 
(llbata et all 11994 [19971 iBellazzini et all 12003ft and in 
M31 ([Perrett et alll200l iMackev et al.ll2010f ). 

For our search, we target NGC 5128 (Centaurus A) 
which is a giant elliptical g alaxy at a distance of 3.8±0.1 
Mpc CH arris et all 12010a!) . While this galaxy is quite 
close, it is still beyond the distance where detecting 
fai nt stellar st r eams is a n easy task . Pre vious work 
bv IMalin etafl (|1983f ) and iPeng et"all ([20021 ) has shown 
there are existing complex shell structures that surrounds 
the central regions out to 15 kpc, attributed to the many 
accretions of s mall neighbouring galaxi es. There are 
also HI shells ( Schimi novich et all TT9 94) and a promi- 
nent warped dust lane ([Graham! 119791 ) is also present 
in the central region where gas and dust is settling into 
the ce ntral potential well. There is evidence for y oung 
stars (jReikuba et all 120011 l2002t lEllis et all 120091) that 
may be aligned with the radio jet. This star formation 
could have been triggered by the jet colliding with cloud 
material brought into NGC 5128 via a merging event. 

NGC 5128 h as 607 confirmed GCs , via radial velocity 
meas ur emen ts (Ivan den Bergh et alJ 119811: iHesser et all 



19841 119861 iHarris et all Il992t |Peng et alj"l2004 
Woodlev etall 120051: jReikuba et all 120071: iBeaslev et ail 
20081: iWoodlev et al ll2010allbl ) and/or resolved images 



from the Advanced Camera for Surveys on the Hub- 
ble Sp ace Telescope ([Harris et all 120061 IMouhcine et all 
l2010bf ). Of these, 563 have measured radial veloci- 
ties. Their mean measurement uncertainty is 42 km s _1 , 
though the range of uncertainties is quite large, with 
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96% of the population having uncertainties < 200 km 
s _1 . The currently knowns GCs are distributed in galac- 
tocentric radius out to 50' (where 1' ~ 1.1 kpc at the 
distance of NGC 5128), although the inner 5' is sparsely 
populated because of the obscuration of the well known 
central dustlane. There are also 780 confirmed PNe wit h 
positional and radial velocity data (|Peng et al.ll2004al ). 
which can be used as stellar tracers in the same way. 

2. GROUP IDENTIFICATION TECHNIQUE 
2.1. Globular Cluster Observational Dataset 

From both the GC and PNe databases, we attempt to 
identify genuine subgroups in the classic way by look- 
ing for objects that are close to each other in both po- 
sition a nd velocity. Our ba sic search approach follows 
that of iPerrett etall <|2003) for the GCs in M31. We 
look first at the GCs. Here the NGC 5128 GC data face 
us with the additional difficulty beyond the situation in 
M31: not only are the positions and velocities seen only 
in projection, but the velocity measurement uncertain- 
ties (averaging more than 40 km s _1 ) are usually larger 
than the expected ~ 20 km s _1 true, physical velocity 
dispersion of a typ ical satellite dwarf in the Local Group 
(|Geha et al.ll2010l) . For this reason, we view our present 
work as only a preliminary step towards finding guen- 
uine physical subgroups. For the same reason, we sug- 
gest that at this stage, only rather simple group-finding 
search procedures are justified. Logical next steps would 
include higher-precision velocity data, and imaging to 
search for faint surface brightness features. 

GC metallicity or color might in principle be addi- 
tional search criteria, as it may be expected that GCs 
within one accreted system may share similar enrich- 
ment processes, but these properties are less useful in 
practice. In our sample of GCs, the metallicities are ei- 
ther from a color transformation or simply base d on a 
color d ivision which we take from iWoodlev et al.l (2007, 
l2010allbh . We have classified the GCs as being either 
metal-poor or blue ([Fe/H] < —1.0) or as metal- rich or 
red ([Fe/H] > —1.0) following previous work on the GC 
system of N GC 5128 (jHarris et al.ll2004bt IWoodlev et al.l 
[20051 120071) or as metal-rich or red if (B - I)> 2.072 
and metal-poor or blue if (B — I)< 2.072, following 
iPeng et al . (2004b). The division between red and blue 
GCs is thus not homogeneously determined and depends 
on a variety of methods as well as a number of different 
photometric studies. Of these GCs, 291 are classified as 
metal-poor, 292 as metal-rich, and 24 GCs have insuffi- 
cient photometry for any of the above transformations. 
In addition, while there is evidence that the metallicities, 
on average, of GCs in dwarf gala xies are more meta l- 
poor than in more massive galaxies (jSharina etalll20inh . 
there is evidence that a range of metallicities of GCs can 
exist within a single dwarf galaxy, for example, in th e 
dwarf galaxies IC10 and UGCA86 (ISharina et al.ll2010h . 
As well, there can exist a spread in the metallic ities of 
the d warf galaxies themselves (see the review of iMateol 
1998). This information, taken together, makes it quite 
difficult to incorporate a defined range of acceptable color 
or metallicity within a sub g roup as a search parameter 
(we note that IPerrett et al.l (|2003| ) also ignored metallic- 
ity as a search parameter). 

We searched for GC subgroups in the outer regions of 



NGC 5128, beyond 10'. Inside this radius, there would 
be a high fraction of false subgroups because the den- 
sity of objects is high, and the biggest source of finding 
false subgroups is due to the accidental 2-dimensional 
projection of objects that are actually far apart along 
the light of sight. Additionally, it is the outer halo that 
should still preserve coherent traces of accreted satellites 
because the much larger dynamical times there prevent 
them from being fully phase-mixed ()Bullock fc Johnston! 
I2005T) . 

As mentioned above, we use essentially a nearest- 
neighbour approach with two linking criteria. For each 
GC outside a 10' galactocentric distance, we calculate the 
separation in projected linear distance between it and all 
other GCs in the catalog (the boundary can therefore ex- 
tend slightly inwards of 10') and group the GCs if these 
are within our imposed linking length. We also included 
the GC velocity and velocity uncertainties in our search 
such that two GCs would match if the range of its velocity 
plus or minus its velocity uncertainty can be subtracted 
from another objects velocity such than their absolute 
difference is less than 20 km s -1 . For example, objects 
with velocities 500 ± 30 km s" 1 and 460 ± 20 km s" 1 
would be linked. To help avoid false group detections 
based on poorly measured velocities, we only considered 
GCs whose velocity uncertainties were less than an im- 
posed limit. We also combined groups that were found 
to have common members, which allows for groups that 
may exist in a chain-like structure. 

A search using a spatial linking length of 1' did not 
yield any subgroups with more than 4 members. When 
we increased the linking length to 2', which is a reason- 
able size for a subgroup from a dwarf galaxy, we did 
indeed find plausible subgroups. With a linking length 
of 2', as well as allowing for velocity uncertainties of up 
to 40 km s -1 , we have detected a plausible subgroup of 
GCs including 6 group members. When we include GCs 
in our search that have measured velocity uncertainties 
up to 50 km s _1 , we increase the number of group mem- 
bers by 1, and call this a possible group of 7 members, 
labelled GC 1. When we increase the allowed velocity 
uncertainty to 80 km s _1 , we end up finding 3 additional 
subgroups, each consisting of 5 members, which we la- 
bel groups GC 2, GC 3, and GC 4. Many of these are 
likely to be included simply because of the large velocity 
allowance and we cannot be sure that these are all real 
groups. If we set our velocity uncertainty limit to 40 km 
s _1 , but increase our linking length to 3', we increase 
the number of members in our subgroup from 6 to 10 
(included as possible members of group GC 1). Many 
of these included GCs may not be real members. The 
possible subgroups are plotted in Figure Q] and are listed 
in Table [T] with columns of the GC ID, group ID, R.A. 
and Decl. in J2000 coordinates, the [Fe/H] value, their 
galactocentric radius, R gc , and their measured radial ve- 
locity, v r , and associated uncertainty, o~ Vr . The radial 
velocities and uncertainties are weighted values deter- 
mined from all previous measure ments in the lit e rature . 
These values are obta ined from IWoodlev et aLl (T 20071. 
IWoodlev et all (l2010ah and IWoodlev et al.l (|2010bir ~ 

The image of NGC 5128 in Figure[T]was obtained with 
the Mosaic //optical CCD camera on the 4 meter Blanco 
telescope at the Cerro Tololo Inter- American Observa- 
tory. It was provided to us in its reduced form graciously 



Subgroups of Globular Clusters 



3 



by Er ic Peng with its reduction described in iPeng et all 
( 2002). We have smoothed the image using a 300 x 300 
pixel box and subtracted it from the original image. To 
enhance the visibility of the low-surface brightness shells 
and arc structures and preserve some of the dominant 
bulge light, we added back 15% of the smoothed image. 
We find that in a general sense, the candidate subgroups 
that we have identified fall in the same radial zones of 
galactocentric distance where the shells are found. How- 
ever, a stronger indication of possible connections is seen 
along the iso photal major axis o f the galaxy (35° and 
215° E of N. lDufour etafl (fl97l ). Toward the south- 
east (lower right in Figure [T]), group GC 1 falls on the 
irregular plume of light, while onthe opposite side to the 
northwest, group GC 3 falls on or close to the two of 
the visible arcs. The remaining two groups, GC 2 and 
GC 4, do not seem to be projected on any known surface 
brightness features. 

2.2. Planetary Nebula Observational Dataset 

We have performed the same subgroup finding al- 
gorithm on the available PNe data in the literature. 
There are 780 PNe with radial v elocity measurements 
(IPeng et al.H2"004al: IHui et alJfl995l) in NGC 5128, which 
have been detected by the [OIII] emission line. While 
there are no velocity uncertainties presented in the liter- 
ature for the PNe, we have assumed a radial velocity un- 
certainty of 20 km s" 1 for each obje ct, which is the typ - 
ical rms velocity error presented in IPeng et al.l ([2004a). 
Using our same subgroup finding criteria as for the GCs 
(linking length of 2', velocity difference of 20 km s _1 plus 
the inclusion of the velocity uncertainties, at a minimum 
distance of 10 ' from the center of the galaxy for the cen- 
tral group member), we have found 4 potential groups. 
These groups are plotted in Figure U and listed in Ta - 
ble [U which provides their ID from IPeng et al.l ()2004afl . 
their group ID, R.A. and Decl. in J2000 coordinates, the 
galactocentric radius, R gc , their radial velocity measure- 
ment, v r , and assumed 20 km s _1 velocity uncertainty. 
Out of the 4 subgroups of PNe discovered, 3 have group 
sizes of 5 — 6, while the remaining group has 9 members. 
These groups have been labelled PNe 1, PNe 2, PNe 3, 
and PNe 4. 

2.3. Properties of Subgroups 

The mean velocity and rms velocity dispersions for 
each group are listed in Table [3l The velocity dispersion 
for each subgroup ranges from 29 — 43 km s _1 for the 
GC subgroups and from 12 — 24 km s _1 for the PNe sub- 
groups, within a reasonable range for a plausible group. 

Our results indicate there are 3 possible overlaps be- 
tween the GC and the PN subgroups. These are PNe 1 
with GC 2 (Group 1), PNe 2 with GC 3 (Group 2), and 
PNe 3 with GC 4 (Group 3). Examining the average 
velocities for each of these subgroups, the average radial 
velocity differences are A(< vqc > — < vpnc >) =39 
km s _1 for Group 1, —136 km s _1 for Group 2, and 
km s _1 for Group 3. The difference in average velocities 
for Group 1 and Group 3 are within the expected uncer- 
tainty range of the measured velocities, which is 42 km 
s _1 for the average GC. The subgroups of Group 2, how- 
ever, do not appear to have strong velocity connections, 
and are likely either not real groups or are overlapping 




Fig. 1. — We show the positions of the potential groups of GCs 
and planetary nebalae in NGC 5128. The group labelled GC 1 
(circles) are the 7 members found with a linking length of 2' with 
allowed velocity uncertainties of 50 km s . We also plot the 3 
additional GC subgroups that are found when GCs with velocity 
uncertainties of up to 80 km s — 1 are included, labelled GC 2, GC 3 
and GC 4 (circles). The 4 plausible PNe subgroups are also shown, 
groups PNe 1, PNe 2, PNe 3 (note one PNe is off the plot), and 
PNe 4 (squares). The locations of the friendless PNe determined 
in Section [4] with N = 10 and crjv = 3 (see text for a detailed 
description of these parameters) are also plotted (triangles). The 
image of NGC 5128 is ~ 25' across with North up and East to the 
left, and objects in the outer regions are not shown (see text). 

by chance in projection. In Sections 12.41 & 13.31 we an- 
alyze the probability that the overlapping can occur by 
chance. 

In very rough terms we can estimate the total luminos- 
ity of the underlying satellite galaxy that would consist 
of either ~ 6 GCs or 6 PNe, typical of what we have 
found in our subgroups. Considering first the GC sub- 
groups, we assume a specific frequency, Sn, which is the 
number of GCs per unit My = —15 galaxy luminosity 
(jHarris fc van den~B crgh 198l[). We adopt S n = 3, typi- 
cal for a normal faint galaxy, though note that IPeng et al.l 
( 2008) show that the range of S n for dwarf galaxies spans 
1-20. With these inputs, we obtain My = —15.8 mag 
and a surface brightness of 24 mag per arcsec 2 assuming 
the galaxy is spherically shaped with a 2' radial extent. 
For a PN subgroup of ~ 6 members, we can estimate the 
bolometric luminosity of the underlying satellite galaxy 
with the luminosity specific PN number, a — Npjy / 'L totl 
where Npn is th e total numbe r of PNe and L tn t is the 
total luminosity (| Jacobvl [i"980l) . ICiardullo et all (|2005f ) 
shows ao.5 ~ 2 x 10~ 9 for galaxies in the range of 
Mb ~ —14 to —16 (please see their Figure 1). How- 
ever, in the case of NGC 5128, the PNe dat a extend at 
least 1.5 mag down the luminosity function ()Peng et al.l 
l2004al) . and we thus estimate ai.5 ~ 11.5 x 10 9 , to ob- 
tain L to t = 5.2 x 10 8 L Q . We obtain My = —16.2 mag af - 
ter a bolometric correction of -0.85 (|Buzzoni et al. 2006). 
Assuming the galaxy is spherical with a radial extent of 
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2', we obtain a surface brightness of 24 mag per arcsec 2 . 
Both estimates of the surface brightness of a potential 
underlying dwarf galaxy are very rough guides due to 
the input assumptions, and we do not think this test is 
capable of ruling out whether or not we should see these 
structures in previously obtained deep images of NGC 
5128. We appear to be finding mild traces of substruc- 
ture that are not hugely above the statistical noise. 

2.4. Scrambling the Velocities 

Our first step towards assessing how real our candidate 
groups in Table [T] are is through a non-parametric test. 
By taking our GC positions in 2-dimensional space, we 
randomly assign each GC a unique and real velocity and 
velocity uncertainty from the GC catalog. We then per- 
form our search for subgroups with a linking length of 
2', a velocity difference of 20 km s~ x plus the inclusion 
of velocity uncertainties up to 40 km s _1 , at a mini- 
mum distance of 10 ' from the center of the galaxy for 
the central group member. Out of 1000 trials, we found 
that 23.2% of cases found one group containing at least 4 
members. In 1.8% of the cases, the trials found 2 groups 
with at least 4 members, and in 1 case (0.1%), there were 
3 groups found. This indicates that it is likely that one 
of the GC candidate subgroups may be a false group. We 
did not find a case, however, in which 4 candidate groups 
arose with this methodology. 

We examined closely the groups in this test that had 
more than one group per trial. In the 18 cases that found 
2 groups, there were 2 trials in which the groups over- 
lapped in projected space. In the trial consisting of 3 
groups, 2 of these groups overlapped. This indicates, 
that while finding more than 1 group in a trial was rare, 
when it did occur, there was an approximate 15% chance 
they would overlap. Since these groups were all from the 
GC population, none of the subgroups that overlapped 
had similar average radial velocities. 

Finally, we randomly assigned each GC a unique and 
real velocity and velocity uncertainty from the GC cat- 
alog, but this time, we incorporated GCs with velocity 
uncertainties up to 80 km s _1 . Our results indicate that 
out of 1000 trials, at least 4 false groups were found in 
1.5% of the cases. This simple test indicates that it is not 
likely that all 4 GC subgroups are false, but it is quite 
probable that some may indeed by false. 

3. SIMULATING GLOBULAR CLUSTER SYSTEMS 

3.1. Searching for False Subgroups 

Here, we assess how often our search procedure would 
turn up false groups of similar apparent quality from 
purely random distributions. It is similarly important 
to examine how often our procedure would successfully 
find a real subgroup within the data. 

To carry out tests of the procedure, we generated 
simulated GC systems in which the particles are lo- 
cated randomly within a spherically symmetric space 
distribution following a 3-dimensional density power law 
p ~ r~ n , with n(3D) = 2.5 for the metal-poor GCs and 
n(3D) = 3.0 for the metal-rich clusters. The particle ve- 
locities were assigned randomly in a Gaussian distribu- 
tion with a dispersion of 150 km s _1 for all GCs, and their 
individual measurement uncertainties were randomly dis- 
tributed uniformly from 5 — 80 km s _1 . All of these pa- 
rameters mimic the NGC 5128 GC system. 



For 1000 realizations, each GC system consisted of 
550 objects, with 50% red and 50% blue subpopulations, 
which were projected from their real 3-dimensional posi- 
tions into 2-dimensions. We then searched for subgroups 
with our finding algorithm, as we have done for the obser- 
vational data in Section |2~T1 Our search criteria allowed 
objects to be group members only if they were within 
a projected length of 2' from each other, have a veloc- 
ity difference of < 20 km s , exist outside of 10' from 
the galaxy's center, and have a velocity uncertainty less 
than 40 km s _1 for each member. Allowing the higher 
velocity uncertainty limit than imposed in most of our 
observational searches of 40 km s _1 will only overesti- 
mate the number of false group members we find in the 
simulations. Out of the 1000 trials, 9 of them detected 
one subgroup. The subgroups consisted of 5-7 members. 
This test therefore suggests that there is a 0.9% chance of 
detecting a false group from our observationally mimiced 
dataset. 

We further this work by generating a GC system that 
has 1300 G Cs, which is th e total estimated population in 
NGC 5128 (|Harrisll2010b[ ). Our simulated system has 650 
blue and 650 red clusters, projected into 2-dimensions. 
In the observed GC system, there are 242 objects be- 
yond our inner search radius of 10' so we randomly se- 
lected 242 objects beyond 10' from our simulated sample 
of 1300 GCs. We repeated our group finding algorithm 
to this simulated system as described above. This simu- 
lation represents the most realistic system we can make 
to match our observed dataset. Out of our 1000 simu- 
lations, there were 3 trials that detected a subgroup of 
GCs, each of these groups consisted of 5 members, sug- 
gesting that we would find a fake subgroup of GCs with 
a 0.3% chance. Both this simulation, as well as our sim- 
ulated group of 550 objects show that the probability of 
finding one fake subgroup of GCs in a system similar to 
our observed dataset is < 1%. None of the trials detected 
more than one subgroup in the simulation. 

3.2. The Addition of Real Subgroups 

It is important to also check how reliable our subgroup 
finding technique is at recovering a real group in our 
observed system. This was tested by inserting a real 
subgroup into one of our simulations of 1300 members 
where no fake subgroup was found. We generated 1000 
real groups by randomly selecting their position in x, y, 
and z in the galaxy, as well as their velocity and velocity 
uncertainty. Each real subgroup had between 4-6 mem- 
bers with a group velocity dispersion of within 20 km 
s _1 . Their individual velocities were selected from a sys- 
tem with 150 km s -1 dispersion. The group was inserted 
between 10-50' from the galaxy's center which ensures 
the real group should be present in the projected system 
in our search regime. The differences in the x and y and 
z positions for the group members was selected within 2' 
of each other. 

Out of the 1000 simulations that we performed, 892 of 
these simulations found all group members and no fake 
interlopers. All real group members were recovered as 
well as the inclusion of fake members in 103 of the sim- 
ulations. In 2 of the simulations, some of the real group 
members were found with fake member interlopers, and 
in 3 of the simulations some of the group members were 
found with no fake interlopers detected. We therefore 
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consider our detection method to be very successful as 
in more than > 89% of cases, we were able to find our 
exact simulation group system and 10.3% of the time, we 
also found all group members however with interlopers. 
When interlopers were present in the detected subgroups, 
96.2% of these occurrances had only 1 interloper, while 2 
interlopers were detected in the remaining cases. These 
results strongly suggest that at least some of our obser- 
vationally detected subgroups are real and they cannot 
all simultaneously be false detections. 

3.3. Testing for the Probability of Random Overlapping 

Groups 

Lastly, we test whether our candidate subgroups might 
be smaller ones overlapped by projection purely by 
chance. To the simulated GC system as a whole, we 
add random subgroups consisting of 4-6 members whose 
members are confined within a 2' region. These sub- 
groups were randomly placed in the galaxy between 10- 
50' before being deprojected into 2-dimensional space. 
Two groups were considered to overlap if any one of their 
members were within a small designated distance from 
any member of another group. The distance chosen did 
not significantly change our findings (considering differ- 
ences between 0.1 — 2'), so we selected 0.5'. 

We performed 7 simulations, each consisting of a differ- 
ent number of subgroups placed in the galaxy, which in- 
creased linearly from 2 subgroups up to 8 subgroups. For 
each of these simulations, we examined, over 1000 trials, 
how many times these subgroups overlapped. The num- 
ber of overlaps ranged from no overlapping subgroups up 
to 5 overlapping subgroups. We did not find any simu- 
lation which had more than 5 overlapping groups. As 
would be expected, the higher the number of subgroups 
present, the higher the number of overlapping subgroups 
are found as a result. These probabilities of having ran- 
dom overlapping groups are listed in Table HJ We find 
that the probability of having two subgroups randomly 
overlap is - 4%. 'if all 8 subgroups in NGC 5128 (4 
GCs and 4 PNe) are real groups, then the probability 
that 3 sets of subgroups randomly overlap is found to 
be ~ 7%. Exluding Group 2, whose subgroups have very 
different average velocities, our results suggest there may 
be 2 sets of overlapping subgroups, Group 1 and Group3. 
The probabilty of 2 sets of subgroups randomly overlap- 
ing with 8 subgroups present is 23%. It is therefore not 
clear if these overlapping subgroups found in Groups 1 
and 3 are real. 

4. SEARCHING FOR STELLAR ORBITS 

Work by iPeng et all (|2002f ) shows a complex but faint 
shell structure threading the inner and mid-halo regions 
of NGC 5128. These shells have been attributed to the 
phase-wrapping of stars from accreted satellites. One 
particular blue elliptical arc was found to be associated 
with a young star clusters (~ 350 Myr), but aside from 
this one example, no stellar streams have been associated 
with GCs or PNe. 

We have begun a search for evidence of stellar streams 
in NGC 5128 by probing the objects within the galaxy 
that may trace the stream. We use the available data for 
the GCs and the PNe, which may be stripped along with 
the stellar stream material during galaxy interactions. If 
these objects are located, they can be used to identify the 



possible orbits of accre ted satellites. To do this, we use 
a friendless algorithm ([Merrett et alj 120031 ) and search 
for evidence of stellar streams using the GCs and PNe 
separately in the analysis below. 

For each object in our sample with a measured radial 
velocity and a velocity uncertainty less than 50 km s _1 
(totalling 411 GCs and 780 PNe), we have found its N 
nearest neighbours. For those N neighbours, we have 
calculated their mean velocity v mean ^N and their velocity 
dispersion, ctjv. We have classified the central object as 
friendless if its velocity is more than n x cr/v from v mea n,N ■ 
Varying the N parameter does not seem to significantly 
alter the results. We do find varying the o~n parameter 
from 3 to 4 changes the number of friendless objects by 
a few, while varying ctjv from 2 to 3 does change the 
number of objects significantly. 

Taking a conservative approach and using o~n — 4, we 
find no friendless GCs when N = 30, 1 when N = 20 
(GC0258), and 1 when N = 10 (GC0578). GC0258 will 
not be consid ered further as it was classified as a re- 
solved GC in Ha rris et alj (|2006|) while a recent radial 
velocity measurement i ndicates this object may be a star 
(jWoodlev et al.ll2010bD . and it is because of this low ra- 
dial velocity measurement that this object is identified. 
GC0578 may indeed be friendless, but it is not possible to 
use one GC to indicate evidence for a stellar stream. Us- 
ing N = 10 and a N = 3 yields 7 friendless GCs (GC0216, 
GC0249, GC0325, GC0506, GC0577, GC0578, GC0580). 

We conduct the same analysis with the PNe dataset us- 
ing <tn = 4 and find no friendless objects when N = 30, 
1 when N = 20 (f05p02), and 6 when N = 10 (4012, 
4285, 5602, f05p02, fl8p28, f42pl0). Using N = 10 
and crjv = 3 yields 20 friendless PNe (4012, 4128, 
4285, 4417, 4511, 5203, 5601, 5602, 5617, 6104, f04pl, 
f05p02, f08 P 50=3002, fl4p016, fl8 P 28, fl8 P 65=1208, 
fl8p67=4023, fl8p83, f42pl0, mosNEpn4). The loca- 
tions of the inner ~ 15 of these are shown in Figure [TJ 

The friendless objects that are located at the largest 
distances from the galaxy's center are likely not indica- 
tors of streams. Rather, there are very few neighbours as 
the incompleteness of planetary nebulae and GCs is quite 
high beyond 20' (see the radial distributi on as a function 
of azim uth for the GCs in Figure 8 of iWoodlev et al.l 
(2010 b]) and the d istribution of the PN in Figure 2 of 
IPeng et alj (|2004al )). This leaves only the inner 20' which 
can be explored for friendless objects, of which there are 
few (again, depending on the parameters chosen in the 
friendless algorithm). We list these objects with the sug- 
gestion that (at least) some of them may trace stellar or- 
bits, and continued searches for stellar streams or orbital 
modelling are necessary to go beyond this preliminary 
work. 

5. CONCLUSIONS 

We have searched the GC and PN systems in the 
nearby giant elliptical galaxy, NGC 5128, for evidence 
of stellar subgroups and stellar streams from merging 
events. Our results indicate there may be up to 4 po- 
tential subgroups of GCs and 4 potential subgroups of 
PNe. Two of the PNe subgroups overlap with two of 
the GC subgroups in position and velocity. In order to 
improve the search for these subgroups in NGC 5128, 
we need higher precision radial velocity measurements to 
more concretely determine if these objects are grouped 
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together. 

By generating GC systems that mimic that in NGC 
5128, we are able to assign the probability of our results 
being a chance encounter due to the projection of these 
objects onto a 2-dimensional plane. The probability of 
finding a fake subgroup within the GC population is < 
1%, while the probablity of having a set of subgroups 
overlap is ~ 4%. These results strongly suggest that at 
least some of these subgroups may indeed be real. 

We have also searched for evidence of stellar streams 
by locating GCs and planetary nebalae that differ from 
its neighbours in velocity space using a friendless search 



algorithm. We present our findings as potential trac- 
ers, however further work in orbital modelling as well as 
searches for very faint surface brightness features would 
be necessary beyond this point to further define the stel- 
lar streams. 



K.A.W. thanks Dr. Harvey Richer and W.E.H. thanks 
NSERC for their financial support. K.A.W and W.E.H. 
thank Dr. Eric Peng for providing the image of NGC 
5128 in Figure [U 
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TABLE 1 

Possible Globular Cluster Subgroups in NGC 5128 





Group ID 


/x.-ri. 


LJCCL . 




tigc 


v r 


rr c 








\ J ZUUU J 


[J zUUU ) 




l'\ 


[km s ) 


[km s ) 


GC0067 




1 


13:24:51.49 


-43:12:11.1 


-1.17 


12.86 


624 


13 


GC0083 


GC 


1 


13:24:56.08 


-43:10:16.4 


-2.46 


10.79 


687 


31 


GC0106 


GC 


1 


13:25:01.83 


-43:09:25.4 


-1.42 


9.52 


688 


16 


GC0114 a 


GC 


1 


13:25:03.37 


-43:11:39.6 


-1.29 


11.41 


605 


46 


GC0117 


GC 


1 


13:25:04.48 


-43:10:48.4 


-0.32 


10.54 


626 


22 


GC0122 b 


GC 


1 


13:25:05.46 


-43:14:02.6 


-1.37 


13.51 


679 


13 


GC0419 


GC 


1 


13:24:55.31 


-43:10:39.3 


-0.35 


11.19 


656 


27 


GC0422 b 


GC 


1 


13:25:03.28 


-43:08:14.4 


-1.65 


8.37 


690 


31 


GC0491 


GC 


1 


13:24:54.98 


-43:11:36.9 


-0.22 


12.05 


670 


19 




GC 


1 


1 ^-9zl- ^fi fil 


/1Q. 1 9-OQ p. 


-1 94 


1 9 


71 ^ 
/ ID 






GC 


1 


13:25:08.94 


-4o:Uo:oo. f 


-1.96 


8.46 


678 


26 


GOU354 


GC 


2 


10 ■ i / ■ no ocr 

13:26:02.25 


-43:08:55.6 


-0.65 


10.03 


457 


38 


GC0598 


GC 


2 


13:26:09.52 


-43:08:52.4 


-0.41 


10.88 


482 


67 




GC 


2 


1Q .' 1 ( ■ -io oo 

13:26:12.82 


-43:09:09.2 


-0.49 


11.50 


545 


60 


/ 1 pnc ,1 /t 

LrOU044 


GC 


2 


13:26:17.27 


/i q . nn .ko n 
-43:09:58.0 


-0.05 


12.65 


A -TO 

473 


39 


GC0393 


GC 


2 


13:26:22.08 


-43:09:10.7 


-1.20 


12.79 


505 


78 


GC0534 


GC 


3 


13:25:55.35 


-42:53:40.2 


-2.30 


9.04 


347 


41 


GC0590 


GC 


3 


13:25:46.96 


-42:52:34.0 


-0.95 


9.28 


355 


66 


GC0532 


GC 


3 


13:25:54.79 


-42:52:57.1 


-0.33 


9.59 


399 


32 


GC0341 


GC 


3 


13:25:58.91 


-42:53:18.9 


-1.34 


9.70 


410 


20 


GC0334 


GC 


3 


13:25:56.59 


-42:51:46.6 


-0.73 


10.77 


403 


16 


GC0518 


GC 


4 


13:25:25.36 


-42:49:27.5 


-1.55 


11.70 


589 


63 


GC0163 


GC 


4 


13:25:15.79 


-42:49:15.1 


-1.44 


12.09 


534 


52 


GC0516 


GC 


4 


13:25:23.31 


-42:48:47.7 


-0.54 


12.38 


526 


11 


GC0220 


GC 


4 


13:25:30.72 


-42:48:13.4 


-1.88 


12.94 


484 


47 


GC0511 


GC 


4 


13:25:18.29 


-42:47:52.4 


-0.65 


13.39 


487 


39 



a Included if velocity uncertainty cut-off increased from 40 km s 1 to 50 km s 
Included member in group if linking length increased from l' to 3'. 



TABLE 2 





Possible Pl 


ANETARY NEB' 


ULAE SUBGROl 


jps in NGC 


5128 




ID 


Group ID 


R.A. 


Deal. 


Rgc 










(J2000) 


(J2000) 


(') (kms- 1 ) 


(km s- 1 ) 3 


ctr513 


PNe 1 


13:26:14.26 


-43:08:22.5 


11.18 


480 


20 


f08pl0=5506 


PNc 1 


13:26:06.19 


-43:09:06.2 


10.63 


462 


20 


5508 


PNe 1 


13:26:13.01 


-43:06:45.0 


10.01 


439 


20 


f07p06 


PNe 1 


13:26:11.33 


-43:09:57.3 


11.89 


421 


20 


f07p08 


PNe 1 


13:26:06.87 


-43:07:36.0 


9.65 


476 


20 


f08p9 


PNe 1 


13:25:56.04 


-43:09:24.7 


9.76 


439 


20 


fl7p34 


PNe 2 


13:26:11.29 


-42:52:39.6 


11.65 


535 


20 


fl8p46 


PNe 2 


13:26:10.83 


-42:54:01.5 


10.64 


527 


20 


fl8p51 


PNe 2 


13:26:07.42 


-42:54:01.4 


10.18 


530 


20 


fl8p55=1217 


PNe 2 


13:26:02.57 


-42:53:07.1 


10.26 


516 


20 


fl8p74=1303 


PNc 2 


13:26:05.35 


-42:51:31.3 


11.84 


522 


20 


5106 


PNc 2 


13:26:11.57 


-42:55:37.7 


9.75 


506 


20 


5108 


PNe 2 


13:26:09.16 


-42:55:20.1 


9.56 


513 


20 


fl8p69=4031 


PNc 2 


13:25:56.78 


-42:53:50.4 


9.05 


499 


20 


fl8p70 


PNe 2 


13:25:58.98 


-42:53:39.0 


9.44 


517 


20 


fl8pl5 


PNe 3 


13:25:41.55 


-42:45:41.5 


15.67 


519 


20 


fl8pl6 


PNe 3 


13:25:38.18 


-42:46:04.4 


15.20 


501 


20 


fl8pl7 


PNe 3 


13:25:39.18 


-42:45:24.2 


15.89 


534 


20 


fl8pl8 


PNe 3 


13:25:34.32 


-42:46:11.5 


15.01 


531 


20 


fl8p40 


PNe 3 


13:25:40.29 


-42:43:44.0 


17.57 


553 


20 


fl8p26 


PNe 3 


13:25:26.79 


-42:47:36.5 


13.54 


507 


20 


ctr624 


PNc 4 


13:24:50.38 


-42:53:22.6 


10.33 


557 


20 


fl9p06 


PNe 4 


13:24:45.76 


-42:51:30.0 


12.31 


588 


20 


fl9p07 


PNe 4 


13:24:44.04 


-42:52:09.9 


12.01 


538 


20 


4608 


PNe 4 


13:24:49.87 


-42:52:46.8 


10.85 


597 


20 


ctr622 


PNc 4 


13:24:53.26 


-42:53:26.8 


9.94 


566 


20 



a For all PNe, we have assumed a radial velocity uncertainty of 20 km s 
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TABLE 3 

Properties of the Subgroups 

Group ID 

(km s- 1 ) (km s- 1 ) 

GC l a 651 25 

GC 2 492 26 

GC 3 382 19 

GC 4 524 34 

PNe 1 453 21 

PNe 2 518 11 

PNe 3 524 17 

PNe 4 569 21 

a These values have been calculated us- 
ing the 7 probably group members of 
GC 1 using a linking length of 2 . 



TABLE 4 
Chance of Overlapping Subgroups 



N Subgroups N Occurances out of 1000 

overlaps 1 overlap 2 overlaps 3 overlaps 4 overlaps 5 overlaps 



2 958 42 

3 869 129 2 

4 783 213 4 

5 631 327 40 2 

6 439 450 99 11 1 

7 358 445 160 36 1 

8 249 441 229 67 13 1 



